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ABSTRACT 

Following the discovery of high-energy (HE; E > lOMeV) and very-high-energy (VHE; 
E > 100 GeV) y-ray emission from the low-frequency-peaked BE Lac (LBL) object AP Li¬ 
brae, its electromagnetic spectrum is studied over 60 octaves in energy. Contemporaneous data 
in radio, optical and UV together with the (non simultaneous) y-ray data are used to construct 
the most precise spectral energy distribution of this source. The data have been found to be 
modeled with difficulties with single zone homogeneous leptonic synchrotron self-Compton 
(SSC) radiative scenarios due to the unprecedented width of the high-energy component when 
compared to the lower-energy component. The two other LBL objects also detected at VHE 
appear to have similar modeling difficulties. Nevertheless, VHE y rays produced in the ex¬ 
tended jet could account for the VHE flux observed by H.E.S.S. 

Key words: gamma rays: observations - Galaxies : active - Galaxies : jets - BE Lacertae 
objects: individual objects: AP Librae 


1 INTRODUCTION 

Blazars are among the most energetic objects in the Universe that 
exhibit non-thermal electromagnetic spectra from radio up to very- 
high-energy (VHE, E>100 GeV) y-rays, with a two-component 
spectral energy distribution (SED) structure in a v/(v) represen¬ 
tation. Multi-wavelength data are of paramount importance to un¬ 
derstand the mechanisms at play in the jet. 

Blazars are divided into two classes: flat spectrum radio 
quasars (FSRQs) and BE Lacertae (BL Lac) objects, the latter be¬ 
ing sub-divided into high-frequency-peaked BL Lac (HBL) and 
low-frequency-peaked BL Lac (LBL). The distinction between 
HBL and LBL classes is b ased on the low-energy peak position 
jPadovani & Giommill995t) . HBL objects present a peak in the UV 
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or X-ray range while the peak of LBL objects is located at lower 
energies (i.e., in optical wavelengths). 

So far, the vast majority of BL Lac objects detected in VHE 
belong to the HBL sub-clasB The SEDs of HBL objects are of¬ 
ten successfully modeled with a synchrotron self-Compton (SSC) 
model, in which the low-energy emission is produced by syn¬ 
chrotron radiation of relativistic electrons, and the high-energy 
component by inverse Compton scattering off the same synchrotron 
photons. HBL are the dominant class of extragalactic objects de¬ 
tected by ground-based Atmospheric Cerenkov Telescopes (ACTs) 
in the TeV y-ray regime. 

Only a few TeV emitters bel ong to the LBL s ub-class and, 
among them, AP Librae (z = 0.049. ljones et alj2009h was recently 


* To keep track of the number of detected object, an up-to-date VHE y-ray 
catalogue can be found in the TeVCat http: //tevcat. uchicago. edu 
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detected by the H.E.S.S. collaboration iAbramowski et al]|2015h 
with a flux of 8.78 ± 1.54stat ± 1.76sys x 10“*^ cm“^s“* above 130 
GeV and a photon index T = 2.65 ± 0.19stat ± 0.20sys matching well 
the spectrum measured by the Fermi Large Area telescope (LAT) in 
the high energy (HE, 100MeV<E<300GeV) range. Remarkably, 
the spectral break between the HE and VHE ranges is the smallest 
ever measured for an LBL object but cannot be expl ained by extra- 
galac tic background light (EBL) attenuation only iSanchez et al.l 
2013t). In this work, VHE and HE data have been extracted from 
Abramowski et alj j2015l) 

After the announceme nt of this detection by the H.E.S.S. col¬ 
laboration lHofmannll20 1 Oh . Swift and RXTE data were taken cre¬ 
ating contemporaneous spectra in X-ray and UV bands. Analysis 
and results are presented in sections Oand ^ Archival ob¬ 
servation by Chandra ( ^2.21 have been analyzed in this work, re¬ 
vealing the first X-ray extended jet for a VHE blazar. At longer 
wavelengths, AP Librae is one of the targets of different monitor¬ 
ing programs such as SMARTS ( 32.3b and the MOJAVE program 
( 32.4t which provide long-term optical and VLBA measurements. 
The VHE detection, together with lower energy-data presented in 
this paper, enable to draw the first complete SED of this source and 
to probe mechanisms at play in LBL objects. The broadband SED 
is then discussed in the framework of different emission models in 

Throughout this paper a ACDM cosmology with Ho = 
71 kms“* Mpc^’, Oa = 0.73 and Dm = 0.27 is assumed, resulting 
in a luminosity distance for AP Libra e of Dl = 21 5 Mpc and a 
linear scale of 0.947 kpc per arcsecond dWri ^l2006t) . 


2 MULTI-WAVELENGTH OBSERVATIONS 

2.1 Swf/i-XRT and RXTE-PCA observations 

X-ray observations of AP Librae during the period of interest were 
retrieved using the HEASARC archive. Four consecutive daily ob¬ 
servations (ObsID 95141) of ^ 3 ks each were car ried out between 
10-14 July 2010 with RXTE i.Iahoda et alJll996h . with a total ex¬ 
posure of ^ 13 ks. The STANDARD2 RXTE-PCA data were ex¬ 
tracted using the ftools in the HEASOFT 6.16 software package 
provided by NASA/GSFC and filtered using the RXTE Guest Ob¬ 
server Facility recommended criteria. Only signals from the top 
layer (XIL and XIR) of Proportional Counter Unit 2 (PCU2) were 
used to extract spectra in the 3 - 50keV range, using the faint- 
background model. The obtained daily light curve has an aver¬ 
age rate of 0.44 counts s“', a variance of 0.03 counts s“^ compati¬ 
ble with its expected variance of 0.02 counts s“^ if the source were 
constant, and a chi-square probability of constancy of 27%, hence 
no variability is present over the span of 4 days. 

During the period of interest, seven obser vations were carried 
out by the Swift mission jBurrows et al.|[2003) . between 20 Febru¬ 
ary 2010 and 16 August 2011 (ObsID 36341005 to 36341011), 
of which one 5ks observation was carried out on 7 July 2010, 
near the RXTE observation. However, the short observation in Ob¬ 
sID 36341009 was skipped. The photon-counting (PC) mode data 
are processed with the standard xrtpipeline tool (HEASOFT 
6.16), with the source and background-extraction regions defined 
as a 20-pixel (4.7 arcsec) and a 40-pixel radius circle respectively, 
the latter being centered nearby the former without overlapping. 
All exposures show a source with a stable average count rate of 
0.12 counts s“'. Also the large 5ks XRT-PC light curve shows 
the source with an average count rate of (0.13 ± 0.02) s“' and 


an r.m.s. of 0.01 s“* for which no variability could be found 
with a 99% c onfidence level upper limit on the fractional variance 
(as defined in I Vaughan et al.ir2003h Evar of 0.95. Using this count 
rate in WebPIMMS from HEASARC, a RXTE-PCA count rate of 
0.6 counts s“* is predicted, compatible with the value actually 
observed of 0.44 counts s“' hinting at the fact that the source was 
probably in the same state during observations of both observato¬ 
ries. Given the low count rate, no pile-up is expected in PC mode, 
which is confirmed by the acceptable fit of a King profile to the PSF 
of all observations. 

Sp ectral fitting o f all ObsIDs was performed with PyXspec 
vl.0.4 {a rnau 3 [T^ . using a response matrix for the combined 
PCA data set generated by the ftool pcarsp vll. 7.1, and dedi¬ 
cated Ancillary Response Functions (ARFs) for all XRT data sets 
generated by xrtmkarf (along with the latest spectral redistri¬ 
bution matrices swxpcOtol2s6_2O110181v014 from CALDB). 
Spectra from all ObsIDs were rebinned to have at least 20 counts 
per bin using grppha, channels 0 to 29 were ignored in the XRT-PC 
data, and only the 3-50 keV range is used in the PCA data. All data 
sets are fit to a power-law model dV/d£ = No(E/E qC^'*^ , where No 
is the normalization factor at a chosen reference energy Eo = 1 keV 
and Tx the photon inde x. Using the Leiden/A rgentine/Bonn (LAB) 
Survey of Galactic HI jKalberla et akllloo^ weighted average hy¬ 
drogen column density of Ah = 8.14 x 10^°cm“^, good fits are 
obtained for the power-law function (P{x^) = 0.18 - 0.9) with a 
photon index of Tx 1.55 on average. All XRT observations were 
also summed, a new exposure file built with ximage, and a new 
ARE for the summed spectrum. This latter spectrum extends up 
to 7keV. Another spectrum was derived this time limited to 1 
count/bin, to allow an extension to higher energies, and was fit¬ 
ted using statistic cstat required in the case of Poisson data. 
The fit parameters are entirely compatible with those obtained us¬ 
ing statistics, but the spectrum extends up to =: lOkeV. All fit 
parameters, along with the unabsorbed 0.3-10 keV flux Eo. 3 _iokcV 
(retrieved for each flux using cflux), are shown in Table[T]and the 
light curve is shown on Figure[T] 

Systematic errors on the Swift-XKT spectra and absol ute flux 
are less than 3% and 10%, respectively l lGodet et al.ll2009 l). while 
PCA-X RT cross-calibration details can be found in Tsuiimoto et al.l 

(l201lh . 


2.2 Chandra observations 

AP Librae was observed by Chandra on July 4, 2003 with a total 
exposure time of 14 ks. The Chandra data reprocessing and reduc¬ 
tion were performed following the standard procedures described in 
the Chandra Interactive Analysis of Observation^ (CIAO) threads, 
using CIAO v4.3 and the Chandra Calibration Database (CALDB) 
version 4.4.6. The data reveal the presence of an extended jet on 
arcsecond scales, which is unique amongst the VHE emitting BL 
Lac class so far. A radio VLA observation was used to align the 
nuclear X-ray emission with the radio core. A registered, exposure- 
corrected and adaptively smoothed image of AP Librae in units of 
phcm^^ s“* px“*, with radio contours overlaid is shown in Figure 
In order to assess to what degree the RXTE and Swift spectra 
need corrections for non-core emission, the spectrum of the jet is 
estimated, with the caveat that this observation is not contempora¬ 
neous with the data set presented here. 
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Table 1. Results of the spectral fitting of all XRT-PC and PCA observations 


ObsID 

Time 

No 

Tx 

P(X^) 

Fo.3-IOkeV 


MJD-5500 

phcm“^ s“^ 


% 

xl0-*3 ergcm-3 s-* 

00036341005 

247.2-247.2 

(9± l)x lO-'* 

1.62 ±0.14 

18 

6.8 ±0.8 

00036341006 

249.5-249.7 

(9 ± 1) X 10-“* 

1.45 ±0.09 

70 

7.6 ±0.6 

00036341007 

384.7-384.8 

(8.3 ± 0.4) X 10-“* 

1.47 ±0.06 

31 

7.2 ±0.4 

00036341008 

608.1-608.2 

(10 ± l)x 10-“* 

1 4Q+0.15 

0.14 

35 

Q Q + l.O 

00036341010 

608.0-608.0 


1.51 ±0.09 

94 

7.8 ±0.6 

00036341011 

609.8-609.9 

(9.3 ± 0.5) X 10-“* 

1 co+U.U/ 

0.06 

60 

7.6 ±0.4 

sum all above 


(9.2 ± 0.2) X 10-“* 

1.52 ±0.02 

99 

7.54 ±0.2 

95141 

387.9-391.9 

L3:';|’;^ X 10-3 

1.74 ±0.16 

91 

5.6 ±0.4 

2009- 
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Figure 1. Light curves of AP Librae in, from top to bottom, VHE, HE, X-rays, LFV, Optical and radio (15 GHz) wavebands. The 4 RXTE observations (ObsID 
95141) were merged together and the seven Swift observations (ObsID 36341005 to 36341011) are shown individually. 
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A spectrum of the jet was taken from a polygon shaped re¬ 
gion which avoids the emission of the core and the ACIS readout 
streak. A core spectrum comes from a 2 arcsecond region centered 
on the core. A background spectrum was extracted from four circu¬ 
lar regions placed to the north and south of the source. The jet and 
background regions are marked in Figure|2] In order to estimate the 
eff ects of pile-up in th e core and jet region, the method described 
by iHarris et al.l ( 1201 ih was used. In the jet region no pile-up was 
found while it was necessary to correct for mild pile-up in the core. 

The spectra of the core and the jet contain =: 4900 and 
200 background subtracted counts, respectively. Both spectra 
were binned to a minimum of 20 counts per bin, and fit in the 
0.5-7.0 keV energy band using an absorbed power-law model in 
XSPEC with the same Ah as in 32.11 The fit of the jet spectrum 
yields a photon index Tjet = 1.59 ± 0.16 and a 2-10 keV unab¬ 
sorbed flux of = (1.07 ± 0.37) X 10“'^ erg cm“^ s“*, with 

^ - 4.4 for 7 dof, or more than an order of magnitude below 

the value measured for the source in 32. II based on the Swift and 
RXTE data, which can hence safely be used as the X-ray flux of the 
core in AP Librae. The jet spectrum is comparable with the spec¬ 
tra of large-scale quasar jets observed by Chandra, which may also 
be sources of relatively i ntense y-ray emis sion (see the discussion 
in ISambruna et al.ll2004l : lKnke et al.ll2008l) . Such a scenario is not 
formally excluded here since an extrapolation of the jet spectrum 
could connect within the experimental errors with either the HE or 
VHE fluxes reported here. Assuming no pile-up, the best power-law 
fit to the core spectrum yields a photon index of Ecore = 1.51 ± 0.03 
and a 2-10 keV unabsorbed flux of = 3.18j;Q[4 x 10“'^ 

erg cm“^ s“'. Using the pileup model in XSPEC, a pile-up cor¬ 
rected spectrum appears however to be softer with Tcore = 
and 2.31 x 10“'^ erg cm“^ s“', with a = 158.4 

for 129 dof. The pile-up model of IPavia (120011) was used in the 
fit of the core spectrum, and the value of the pile-up parameter 
a > Q indicates that the fit is indeed affected by this. However, 
it was not possible to obtain an error estimate on a, and hence we 
also do not have an error estimate on the unabsorbed and pile-up 
corrected flux. Due to pile-up effects, the fit results for the core 
should be treated with caut i on. Th is extended X-ray jet was first re¬ 
ported bv iKaufmann et all (l2013h . Our results differ slightly, prob¬ 
ably because we used different extraction and background regions, 
and Kaufmann et al. did not take into account the above-mentioned 
ACIS readout streak. 


2.3 Swi/f-UVOT and SMARTS observations 


All of the available archival data taken on AP Librae with the ultra¬ 
violet and optical telescope (UVOT) on the Swift satellite were ana¬ 
lyzed. This comprised 35 exposures taken between April 2007 and 
July 2010, 13 of which occurred during the time frame with which 
this paper is concerned (see Eigure[T). After extracting the source 
counts from an aperture of 5.0"radius around AP Librae and the 
background counts from four neighboring regions, each of the same 
size, the magnitu des were compu t ed usi ng the uvotsource tool with 
calibrations from iBreeveld et ^ (|201 ll ) . Thes e were converted to 
fluxes using the values from Poole et ^ ( 1200 8h after correction for 
extinc tion following the procedure a nd R„ value of iRoming et al.l 
(l2009l) . The values of a and b from [Roming et al.l |2009l) . com- 

' m 


puted following the procedure ofICardelli et al 


1 19891), were used. 


The E(B - V) value from ISchlaflv & Einkbeinej ( 1201 ih , accessed 
through the NASA/TPAC Extragalactic Database, was used. Results 
are summarized in Tablej^ 



Figure 2. Adaptively smoothed, exposure corrected X-ray image obtained 
by Chandra in the 0.5-7 keV energy band in units of ph cm“^ s“* px“* with 
the pixel size of 0.492", where one arcsecond corresponds to 0.947 kpc on 
linear scale. Overlaid are 1.4 GHz radio emission contours at 2" resolution 
from reprocessing archival VLA data (program AB700; C.C. Cheung, pri¬ 
vate communication). The flux densities of radio contours increase by a fac¬ 
tor of 1.5, starling from a value of 5 times the rms noise equal to 1.82 X10“^ 
Jy beam^L The small red arrow in the radio core shows the orientation of 
the milliarcsecond scale radio jet seen in VLBA (see e.g. iLister et al.l 19981 : 
IZensus et al .I 2 OO 2 I) . and references therein, for a discussion on the radio jet 
at different scales). A black dashed polygon delimits the region used to cal¬ 
culate the jet spectrum, the dashed circles are the background regions and a 
black dashed line indicates a location of an ACIS readout streak. 


AP Librae was ob served in c o ntext of the Yale 
Fermf/SMARTS projecjl dBonning et al.l l2012h . Observations 
were performed in the B, R, J and K bands between Eebruary 
27, 2011 (MJD 55619) and March 3st, 2013 (MID 56739) and 


shown on Figure [T] The number of observations and the mean 
magnitudes are given in Table |2] together with the corresponding 
fluxes. Magnitude s have been corrected for G alactic absorption 
using values from ISchlaflv & Finkbeinei ] (|201lb and converted in 
flux units using the Bessell zero points ( Besselll[T990h . 

The host galaxy of AP Librae is bright and therefore the con¬ 
tribution from starlight must be taken into account to estimate the 
non-thermal flux from the core in the near-infrared to UV band. 
The dereddened near-infrare d and optical measur ements of AP Li¬ 
brae reported in Figure 1 of iFalomo et al.l (Il993h . where the total 
emission was modeled with a giant elliptical galaxy template and a 
superposed non-thermal power-law continuum, are given for illus¬ 
tration in the composite SED of Eigure The synchrotron emis¬ 
sion probably peaks in the optical- to near-IR range, since the spec- 
tral index for AP Lib rae in that range is otirop = 0.95 ± 0.10. In 
iHvvdnen et al.l (l2007h . the fluxes in the B and U bands were calcu¬ 
lated for the host galaxy and the core. The fractional contribution 
of the latter was =; 42% in the B band and 69% in the U band. 
At higher energies the emission from the core accounts for an even 
higher perc entage. To take th is result into account, the host galaxy 
template of ISilva et alj d 19981) has been used and with a normaliza¬ 
tion ajusted to fit the data. 


2.4 MOJAVE 

The parsec-scale structure of the radio jet of AP Librae has been 
monitored throughout the past decade as part of the MOJAVE pro- 
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Table 2. Summary of the Swift-WCfl and SMARTS results. Columns 1 and 3 give the 
filter and corresponding energies and the second column gives the number of observations. 
Magnitudes (column 4) are not coiTected for Galactic absorption. The last column gives 
the connected flux. 


Filter 

-Vobs 

Energy 

[eV] 

magnitude 

Elux 

[10“'* erg cm“^ s“*] 

SMARTS: 

K 

25 

0.56 

11.63±0.37 

1.96±0.78 

J 

79 

0.99 

12.76±0.17 

3.26±0.54 

R 

74 

1.77 

14.53±0.18 

1.88±0.34 

B 

29 

2.86 

15.85±0.13 

1.35±0.18 

Swift-VYOT. 

V 

3 

2.30 

15.18 ±0.04 

2.43 ±0.10 

B 

3 

2.86 

15.94 ± 0.04 

1.69 ±0.06 

U 

7 

3.54 

15.68 ± 0.04 

0.96 ± 0.03 

UVWl 

6 

4.72 

15.88 ± 0.05 

0.63 ± 0.03 

UVW2 

10 

6.12 

16.12 ±0.05 

0.57 ± 0.02 

UVM2 

3 

5.57 

16.09 ± 0.06 

0.55 ± 0.03 


grarr|3 (Monitoring of Jets in Active galactic nuclei with VLBA 
Experiments) with the Very Long Baseline Array (VLBA) at a fre¬ 
quency of 15 GHz. The VLBA data have bee n calibrated and an - 
alyzed following the procedures described by iLister et^ bOOgI) . 
The source shows a bright, continuous inner jet region with a bright 
jet core, i.e apparent jet base, extending towards the South. At a res¬ 
olution of typically ^ (1.5 x 0.5) milli-arcsecond (mas), the core is 
not clearly separated from the inner jet. Elliptical Gaussian compo¬ 
nents were used to model the brightness distribution and to deter¬ 
mine radio flux densities of different emission regions within the 
source. Eor the comparison with higher-energy multiwavelength 
data, we focused on the inner 1.5 mas (=: 1.41 pc) region, which 
could typically be modeled with 2-3 Gaussian model components. 
We have used different models with circular and elliptical model 
components and tested the formal statistical model-fitting uncer¬ 
tainties of the total flux density, which turn out to be much smaller 
(< (1 - 3) %) than the absolute calibration uncertainty, which can 
be conservatively estimated to be of the order of < 10 %. 

The 16 MOJAVE observations from MJD 53853 to 55718 do 
not show sign of significante variability in the VLBI core region. 
Figure |3] shows the value of 1.48 Jy of the radio flux density, aver¬ 
aged over the full observations, from the inner 1.5 mas jet core. 


3 DISCUSSION 

3.1 The radiative components 

The composite SED of AP Librae is shown in Figure [3] To¬ 
gether with the MOJAVE, SMARTS, Chandra, Swift-VWOT, Swift- 
XRT, RXTE, Fermi-LAT and H.E.S.S. data analyzed in this work, 
archival data from NED are reported. In the 30-353 GHz band, the 
Planck measu rements from the Early Release Co mpact Source Cat¬ 
alog (ERCSC. IPlanck Collaboration et al.l201lh are in good agree¬ 
ment with the ar chival data as are th e Wide-field Infrared Survey 
Explorer (WISE, IWright et al.l20I(]|) data in the bands 3.4, 4.6, 12, 
and 22 yum. 


http://www.physics.purdue.edu/astro/MOJAVE 


An extrapolation of the hard X-ray to the optical-UV power- 
law spectrum reported here underestimates the simultaneous 
UVOT flux by at least 2 orders of magnitude, though the steeply 
falling UV spectrum possibly connects with the onset of the XRT 
spectrum. This indicates the presence of an inflection point in 
the SED widely attributed to a transition from synchrotron to IC 
dominated radiation. This feature shows that the Compton com¬ 
ponent of AP Librae is the broadest ever observed in any blazar, 
spanning more than 10 decades in energy from O.lkeV to 
1 TeV. Indeed, only two other objects of the same class as AP 
Librae, and hence with broad Compton components, have been de¬ 
tected at VHE energies so far: BL Lacertae (z = 0.0 69), the first 
LBL object to be proved as being a VHE emitter jAlbert et al.l 
l2007h. and S5 0716-H7 14 (z = 0.310) following an optical trigger 
dAnderhub et alJ2009l) . The observed VHE spectrum of the former 
is not as energetic as AP Librae, and the X-ray spectrum of the 
latter appears to still belong to the synchrotron component. 

An empirical characterization of the two radiative compo¬ 
nents, through a third degr ee polynomial fit of each hump in vFy 
representation (as in, e.g.. lAbdo et alJEoiOl) . is used to estimate 
the synchrotron and IC peak energies. The values of the param¬ 
eters obtained from a fi0 are given in Table and the results 
are represented in the composite SED of Figure [3] As mentioned 
above, the SMARTS and the Swift-WGT measurements in the V, 
B and U were no t used in the fit of th e synchrotron peak as well 
as the data from Ipalomo et al.l l ll993h . The position of the syn¬ 
chrotron peak is then estimated to be Fs.pcak - 0.18 ± 0.06 eV, 
which is compatib le with the value of Fs.peak = 0.26 eV derived by 
lAbdo et ^ J2010h on a different data set. The same authors esti¬ 
mated Fic.pcak = 2.6 /;j'^ GeV for AP Librae in Table 13 based on 
a strong correlation of Fic.peak with the HE photon index The, as 
expressed in their Equati on ^ Using the photon index found by 
lAbramowski et al] J2015li. which is based on an order of magni¬ 
tude larger data set, yields a lower but still compatible value of 


^ The EBL absorption has been taken into account in the fit. 

® The quoted uncertainty, not given in their table, is derived from their esti¬ 
mation of an error of 0.7 associated with the estimation of the log of Fj;,, peak 
in Equation 5. 
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£^ic,pcak = 0.9^,!, j GeV and this value was constrained to be be - 
low 1 GeV by fitting the HE-VHE data dAbramowski et al.ll2015h . 
The polynomial fit presented here yields a much lower value of 
£^ic,peak = n^g'^MeV, which can be attributed to use of the entire 
SED. This is the lowest IC component peak ever measured for a 
TeV-emitting blazar. 

The third degree polynomial also provides a straightforward 
estimation of the curvatures and kiq at the peak positions ifs.peak 
and Eir naak. respec tively, which pertain to the widths of each hump. 
IPaggi et ^ d2009ll show that a relation = 2kic is expected in a 
pure Thomson scattering regime, using a logparabolic parametriza- 
tion of each of the two humps generated by a single zone homo¬ 
geneous SSC model, while Ks = kic/5 in the Klein-Nishina (KN) 
regime. The curvatures found here for AP Librae yield a surpris¬ 
ing Ks - 6.6/fic, emphasizing the broadness of the IC component 
compared to the synchrotron hump, which is hardly possible to re¬ 
produce with simple radiative models. 


3.2 Radiative scenarios 


In a one zone homogenous SSC framework, electrons produce syn¬ 
chrotron photons which are upscattered through the IC mechanism 
by the same electrons to generate the HE and VHE photons. If this 
upscattering occurs predominantly in the Thomson regime up to 
the peak energy, then it becomes possible to constrain the product 
of the magnetic field B and the Dop pler factor 6 for a sing le zone 
homogenous SSC model (following iTavecchio et alJI 19981 , Equa¬ 
tion 4): 


Bd = (l+z)- 


8.6 X 10’£2 


s, peak 


^ic, peak 


( 1 ) 


where the peak energies are expressed in eV. Using the range for 
£s,peak found previously and Eic.peak = 17 MeV yields B6 = 0.17 G. 
The value of the break Lorentz factor yj of the underlying electron 
distribution can also be derived from the ratio of the peak emission 

energies as - 8-5 x 10’. 

Assuming now that the observed synchrotron radiation does 
not exceed 0.1 keV, (i.e. the lowest energy bin in the XRT spec¬ 
trum), which is more likely to belong to the onset of the IC com¬ 
ponent, then this constrains the maximal Lorentz factor of the 
underlying electron population through the maximum synchrotron 
energy 


E ~ A 

^s,max — /m 


BSm^A 

'Scr(l+Z) 


< O.lkeV, 


where B^ = 4.414 x 10'’G is the critical magnetic field leading to 


r^ax < lo’B-'/’r'/’. 


( 2 ) 


Using Equationl^ and Equation.[T] then yields < 2.4 x lO’, 
which is consistent with being a factor -y/Es, max/Tfs, peak higher than 
yt as expected. 

Supposing that electrons with an apparent energy of dy^ax 
have sufficient energy to upscatter photons to at least the maximal 
observed Compton energy £ic,max = 1 TeV, then the Doppler factor 
is constrained to a reasonable value of d > 10. If the scattering of 
0.1 keV photons occurs in the Thomson regime, the Doppler fac¬ 
tor should be such that 4y„,ax x 0.1 keV < 6m,, A. Using the value 
for ymax found above leads to an unusually high value of d > 163. 
If however the scattering occurs in the KN regime for these high- 

r 2 

est energy seed photons, then ific.max = "f^Tmax which, combined 
with the above constraint (Equation|^ on y^ax. then yields 


B6-^ < 2.3 X 10^’ G, (3) 

from which follows, using the above constraint 6 > 10, a reasonable 
constraint of 6 < 2.3 x 10“’ G. In Appendix I aI similar conclusions 
are drawn for an arbitra ry type of seed photons. 

iLister et ^ ( |20 1 31) measured an apparent superluminal motion 
6.4c. This is compatible with d > 10 for a viewing angle below < 5 
degrees and with 6-20 for 1.7 degrees. 

Going further by assuming that photons with energies up to 
Tfic.peak are produced in the Thomson regime, and the 1 TeV pho¬ 
tons in the KN regime, then Equation[T]and Equationj^can be com¬ 
bined to give 6 < 2 X 10“’ G regardless of the value of d. 


3.3 Application of an SSC model to the SED 

The time-averaged SED of AP Librae was modeled with a cano ni- 
cal one zone homogeneous SSC model l lSand & Grindlavll985h . A 
spherical region of size R, with an electron distribution Ndy), mov¬ 
ing with a bul k Doppler factor d, is fil led uniformly with a magnetic 
field B. As in iTavecchio et al.l ( l2010h . Ne{y) is described by a bro¬ 
ken power-law of index S i between y = 1 and yb and S 2 between 
yb and y^ax- The electrons lose their energy by synchrotron emis¬ 
sion, producing a field of photons which become the targets for the 
same electron population through the IC pro cess. The K N effects 
are taken into account using the Jones kernel ( |joneslll9^8l) to com¬ 
pute the IC cross-section. 

A tentative model is shown in Eigure[3] where the shape of the 
electron distribution (5i, ^2 and yb) is constrained by the observed 
synchrotron component. The remaining parameters {R, B, 6, and the 
total number of electrons Ve,tot) are adjusted to reproduce the onset 
of the Compton component in the X-rays. The obtained parame¬ 
ters and model curves, as given in Table |4] and Eigure[3] respec- 
tively (together with t he model parameters and curves derived by 
ITavecchio et al.l ( 1201 (J) for comparison) obey the constraints found 
in 33.21 Not surprisingly, the broad IC component of the SED is 
difficult to reconcile with the synchrotron distribution using such a 
simple model, for which strong indications were already presented 
in 30 

The SSC calculation reproduces well the lower energy part of 
the SED, up to the X-rays, but the spectral prediction in the Fermi- 
LAT energy range is much softer, as well as about one order of 
magnitude above the observed HE flux. The direct consequence 
of the broadness of the IC component is that the H.E.S.S. flux is 
largely underestimated. Directly linked to the electron distribution 
and to the well measured synchrotron component, this shape can 
only be affected by the KN effects, which tend to soften the spec¬ 
trum leading inevitably to even larger disagreements. 


3.4 VHE y Rays from the extended jet? 

As seen in the previous sections, one-zone SSC models can- 
not reproduce the broadband S ED of AP Librae. However, 
[Sottcher. Dermer, & Einkel (I2OO8I) proposed that the Compton¬ 
scattering of the cosmic microwave background (CMB) by elec¬ 
trons in an extended kpc-scale jet could make VHE y-rays. This 
model was suggested to explain the har d VHE spectrum from 
lES 1101-232 as observed by H.E.S.S. dAharonim et al.l[20061 
I2OO7I) . when EBL attenuation was taken into account with the mod¬ 
els available at the time. AP Librae has an extended kpc-scale jet 
resolved in radio (see Pigures[3and|4ll and X-rays (see section lT^ . 
and it has long been thought that the Compton-scattering of CMB 
photons could produce the X-rays observed from these extended 
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Table 3. Parameters of third degree polynomial function de¬ 
scribing the low and high-energy component of AP Librae. 
The function is of the form f{E) = po + Pi logio(^/®^) + 
Pi log?o(£/eV) + p3 log^gCfi/eV) 


Energy Range (eV) 

PO 

Pi 

Pi 

Pi 

3 X IQ-'^ - 50 

-10.79 

-0.52 

-0.41 

-0.048 

50- 10'3 

-13.36 

0.82 

-0.068 

-0.001 


Table 4. Parameters of the SSC model presented in this work and from 
IXavecchio et alJ j201Qh . For both models, ymin = 1 was used. 


Model 

Lb 

Tmax 

Si 

S 2 

^e, tot 

B 

R 

S 


10^ 

10“* 



10” 

10-2[G] 

10*** [cm] 


This work 

1.1 

2.3 

2 

4.9 

5.4 

0.9 

3.5 

20 

Tavecchio et al 

2.0 

5 

2 

4.9 

0.4 

1.2 

1 

40 


v[Hz] 

10® 10" 10^^ 10^® 10’^ 10’® 10^^ 10“ 10^® io^1o^® 



Figure 3. The broadband SED of the LBL AP Librae. The orange triangles 
come from the Planck ERCSC. The orange diamonds are WISE measure¬ 
ments. Blue points are, from low energy to high energy, MOJAVE(15GHz) 
, 5w//?-UVOT (2.30-5.57 eV), SMARTS (0.56-2.86 eV), Swift-XKT/RXTE 
(3-50 keV). Grey points and butterflies are Fermi-LAY for the quiet (cir¬ 
cle) and flare (squa re) periods (0.3-30 0 GeV) and H.E.S.S. (E>100GeV) 
measurement from iAbramowski et 1i]l2015l) . Swi^UVOT, SMARTS data 
are corrected for Galactic extinction and X-ray data are corrected for Ah ab- 
sorpt ion. Light gray data ^ e taken from NED. The dai'k gray triangles come 
from lpalomo et al] Il993b . The red point is the ra dio flux of the ext ended 
jet. The orange line is the host galaxy template of ISilva et al.l il998b . The 
fit with two third degree polynomial functions, not con'ected for EBL, are 
shown with a green line (see ED- The red butterfly is th e Chandra spec¬ 
trum f rom the jet. The dashed line is the SSC model from lTavecchio et~51 
iioID) whereas the red line is the model obtained in this work (see Table 
0 . 


in iFinke et all (l2008h . The compact component can explain the ra¬ 
dio, optical (not including emission that is clearly from the host 
galaxy), X-ray, and the lower-energy Fermi-LAT y-ray data. The 
extended component can explain the extended radio and X-ray 
data, as well as the highest y-ray emission detected hy the LAT 
and H.E.S.S. A douhle-hroken power-law was used to describe the 
electron distribution in the compact component, while only a sin¬ 
gle broken power-law was needed for the electron distribution in 
the extended component. Parameters in the compact component 
are broadly comparable to synchrotron/SSC modeling results for 
other BL Lac objects, including the jet power in electrons being 
sever al orders of magni t ude greater than th at in t he magnetic field 
(e.g., Finke et al.ll2008l : lAbdo et ^l2011allDldll3 : lAliu et aP 1201 31 
l2014j bl). The extended jet is much closer to equipartition between 
electron and magnetic field density by design; a model out of 
equipartition would still be able to reproduce the data. These pa¬ 
rameters are also close to previous results for modeling extended 
jets, although the ma gnetic field is a bit low er than usual (typically 
found > 1 fiG', e.g., Tavecchio et al.ll2007h . This may be because 
previous ICCMB models of extended jets are for FSRQs, rather 
than BL Lac objects. One hypothesis can be that the magnetic fields 
in extended jets of BL Lac objects are lower than those in the ex¬ 
tended jets of FSRQs. 

It should be noted that the ICCMB model for explaining the 
X-ray emission from extended jets is controversial. It could be that 
X-rays are instead produced by synchrotron emi ssion from another 
population of electrons in the extended jet fe.g.. lAtovan & Dermej 
I 2 OO 4 I ; iHardcastQ l2006h . In this altern ative framework, HE and 
VHE emission is unlikely. Recently, iMever & Georganopoulod 
(I2OI4I) used Ferrai-LAT observations to rule out the ICCMB model 
for the X-ray emission from the extended jet in the FSRQ 3C 273. 


jets (e.g.. iTavecchio et al.ll200d ; ICelotti et alj|200lh . Therefore, it 
seems natural to apply this model to AP Librae, to see if the ex¬ 
tended jet emission could plausibly make up the VHE y rays. Thus, 
the broadband SED of AP Librae has been modeled with a com¬ 
pact, synchrotron/SSC model based on iFinke et al.l ( l2008h . and an 
additional component from the extended jet, emitting synchrotron 
and inverse Compton-scattering of CMB photons (hereafter IC¬ 
CMB). 

The result of this model is shown in Figure [4] with the model 
parameters in Table |5] The model parameters are fully described 


3.5 Comparison with other LBL objects 

The SEDs of LBL objects detected in VHE y-rays challenge single 
zone homogeneous SSC radiative models, which usually reproduce 
reasonably well the time-averaged SEDs of the HBL class. 

The most com plete simultaneous c overage of the BL Lacertae 
was established bv lAbdo et al.l ( 1201 Idh during a multi-wavelength 
campaign including the Fermi-GKT and the X-ray observatories 
mentioned in this study for the high-energy part. The X-ray spec¬ 
trum during that campaign was soft, indicating that its origin was 
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Table 5. Model parameters for the SED shown in Fig.|4] The redshift z is 0.049. 


Parameter 

Symbol 

Compact component 

Extended Jet 

Bulk Lorentz Factor 

r 

20 

8 

Doppler factor 

Sd 

20 

8 

Magnetic Field IG] 

B 

0.05 

5.6 X 10^’ 

Variability Timescale Is] 

L 

3.0 X 10“' 

1.35 X 10" 

Comoving radius of blob [cm] 

K 

1.7 X 10'*’ 

3.08 X 10” 

First Electron Spectral Index 

Pi 

2.0 

2.0 

Second Electron Spectral Index 

P2 

3.0 

4.0 

Third Electron Spectral Index 

Pi 

4.2 


Minimum Electron Lorentz Factor 

Xjifn 

1.0 

2.0 

Break Electron Lorentz Factor 1 

y'brKi 

2.8 X 10^ 

4.9 X 10" 

Break Electron Lorentz Factor 2 

ybrk,l 

6.8 X 10^ 


Maximum Electron Lorentz Factor 

y max 

1.0 X 10’ 

2.0 X 10^ 

Jet Power in Magnetic Field [erg s“^] 

Pj.B 

2.2 X 10''’ 

1.4 X 10"" 

Jet Power in Electrons [erg s“^] 

Pie 

1.7 X 10"^ 

2.8 X 10"" 


v[Hz] 

10® 10” 10" 10" 10” 10" 10®' 10®® 10®= 10®l0®“ 



Figure 4. Same as Fig.[^ The red line is the results of the SSC model from 
the compact component and the blue line is the flux originating from the 
extended jet; parameters are given in Table|^ Purple line is the sum of both. 


synchrotron radiation rather than Comptonized photons, making 
for a wider synchrotron vF, distribution than is reported here for 
AP Librae. The difficulty in this case for modeling BL Lacer- 
tae was that the simulated SED required the energy densities to 
be f ar from equipartition . However, a 1997 Beppo-SAX observa¬ 
tion jRavasio et al.ll2003l of BL Lacertae showed a clear IC origin 
for the X-ray radiation, yielding a narrower synchrotron distribu¬ 
tion, for which the SSC model failed to reproduce a reasonable 
(non-simultaneous) HE spectrum, and an external contribution was 
added. 

The broad Compton distribution of S5 0716-t-714, with emis¬ 
sion up to 700 GeV, is either an order of magnitude below the 
best SSC model prediction from lAnderhub et"^ ^20091) . or is too 
wide if the Fermi-hJAI spectru m constrains the flux at Eic.peak (Fig¬ 
ure 6 in iTavecchio et al.|l201(]| ; see also the similar situation for BL 
Lacertae in the same Figure). Note that the HE and VHE data were 
not taken simultaneously in these two LBL objects. 


4 CONCLUSIONS 

Contemporaneous observations of AP Librae with many currently 
available space- and ground-based instruments have been pre¬ 
sented. The data have revealed the broadest Compton distribution 
of any known blazar to date, which spans from X-ray to TeV ener¬ 
gies. 

The SED of AP Librae is difficult to reproduce with a sin¬ 
gle zone SSC model: the steep UV spectrum, probably synchrotron 
emission, does not connect smoothly with the X-ray spectrum, 
which is underestimated by an order of magnitude if a match is 
re quired with the HE y -ray spectrum (as was also pointed out 
bv lTavecchio et al1l2010h . If a match is required with the X-rays, 
the Fermi-hfiH spectrum is then largely overestimated. The new 
H.E.S.S. spectrum further complicates the situation, as none of the 
previous constraints allows this SSC model to reach the VHE do¬ 
main, even assuming a predominantly Thomson scattering regime 
which yields Compton components roughly twice as large in vFy 
as the synchrotron component. There are ways out of the conun¬ 
drum but at the cost of increased model c omplexity. An exam ple 
is blob-in-jet model, recently proposed by iHervet et al.l ( 1201 5h to 
reproduce the SED of Ap Librae. Another possiblility is a model 
where electrons also upscatter soft photons originating outside of 
the high-energy emission site. It has been shown in this work that 
VHE 7 rays from the extended jet, seen in X-ray, can be produced 
and can explain the H.E.S.S. spectrum. 

AP Librae is the third of VHE detected LBL-type object for 
which single zone SSC models fail to reproduce the SED, and is 
currently the only BL Lac type object combining VHE emission 
and a resolved X-ray jet. The LBL class of VHE emitting objects 
proves to be an interesting laboratory to test radiative model sce¬ 
narios, and perhaps to identify parameters on which the LBL-HBL 
sequence could depend. 
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APPENDIX A: CONSTRAINTS FOR AN ARBITRARY 
FIELD OF SEED PHOTONS 


In leptonic class models, the inverse Compton process is respon¬ 
sible for the high energy part of the SED. The seed photons origi¬ 
nate either from synchrotron radiation produced within the jet (SSC 
models) or from a source outside of the jet (external Compton mod¬ 
els). In the latter case, the sources can be either the broad-line re - 
gions or the dust torus dSikora et al.ll994l : lBlazeiowski et alj2000li . 

The peak observed energy £s of an electron with Lorentz fac¬ 
tor 7 is given by 


EJm^c^ 


6y^B 
(1-1- z)B„ 


and the Compton-scattered photon energy by 


EidnieC- = 


V<eed 

(1+z)' 


where the energ}0 of the seed photons is e^ed (respectively = 
dfsecd in the jet’s frame). 

Efficient Compton scattering will occur only for electrons be¬ 
low the KN limit: 


7 < (4<eed)^‘- 


(Al) 


This KN limit means that Compton-scattered photons will be 
mainly restricted to energies: 


EicjlfleC ^ 


16(1 -Hz)e' 


The synchrotron photons produced by the electrons having the en¬ 
ergy (4e'|.j.j)“* have a peak energy given by: 


EsIm^A 


SB 

16(1 +^)<eedfi-' 


Combining the last two equations with the constraints on maximal 
values for E^ w 0.1 keV and E^c ~ 1 TeV derived from the obser¬ 
vations yields: 


6 B 
70 “ 10-2G’ 


(A2) 


which requires either an unusually high Doppler factor, or an un¬ 
usually low magnetic field. If the 1 TeV photons are produced by 
IC scattering in the KN regime, Eg. lAIl becomes 


7 > (4<eed)"‘ 

and the observed photon energy is dTavecchio et al.|[l998l) 

Sy 


E^lmi,d = 


(1+z)’ 


Then Ea. lA2l reads 


6 B 
17.5 “ 10-2G 


(A3) 


which is a reasonable constraint. Note that this calculation applies 
no matter what the seed photon source is (broad-line region or dust 
torus or synchrotron photons produced within the jet), illustrating 
the difficulties of either radiative scenarios to account for the main 
SED features of AR Librae in the Thomson regime. 


APPENDIX B: CANDIDATES FOR VHE OBSERVATIONS 

The detection of AR Librae by the H.E.S.S. telescopes has revealed 
the broadest IC component for a blazar with a peak position at very 
low energy. Unfortunately, only a handful of LBL-type objects have 
yet been detected at VHEs. To decide if AR Librae is a special case 
or a typical representative of the LBL class, other LBL objects have 
to be observed by Cerenkov telescope and detected at VHE. 

Due to their limited field of view (»: 5 deg), an extra-galactic 
survey performed by Cerenkov telescopes is not possible yet. As 
a consequence, good targets for observations have to be found 
based on multi-wavelengths data. In this appendix, six LBL-type 
objects, present i n the second catalog of Fermi sources (2EGL, 
iNolan et al.l20Tll were selected based on their possible VHE emis¬ 
sion. The 2EGL best fit power-law, measured in the 100 MeV- 
100 GeV band, was extrapolat ed above 200 GeV and E BL cor¬ 
rection was made based on the lEranceschini et al] (l2008l) model. 
The redshift informati on was extracted either from the second cat¬ 
alog o f AGN (2LAC. IXckermann et ahlEoi ll) or from ishaw et al.l 
(l2013h . Sources without redshift measurement were excluded and 
only sources classified as a BL Lac of the LBL class were retained. 
Note that AR Librae appeared to be the first on this list when build¬ 
ing it. 

The names of six candidates, ranked by predicted flux above 
200 GeV, are given in Table IB 11 Eor illustration, their SEDs, 
built from archival data using the ASDC SED buildej^ are pre¬ 
sented in Eigure Im1 Two out of the six sources can be ob¬ 
served by H.E.S.S. and five by the northern facilities (VERI¬ 
TAS and MAGIC). Despite its location and with a redshift of 
z = 0.424, the source 2EGL J0738.0-1-1742 can be well suited for 
H.E.S.S. II telescope observations given the lower energy threshold 
(50 GeV) of the instrument. The redshifts of 2EGL J1150.1-I-2419 
and 2EGL J1150.1-I-2419, found in the 2LAC, were not con¬ 
firmed by IShaw et ^ ( l2013h . Eive out of 6 are also present in 
the first Eermi-LAT Catalog o f Sources Above 10 GeV (lEHL, 
The Eermi-LAT Collaborationll2013l . see Eigure ISTt . 


^ The notation E = emgd is adopted here. 


http ://tools. asdc. asi. it/SED/ 































From Radio to TeV: The surprising Spectral Energy Distribution ofAP Librae 11 
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Figure Bl. SEDs for the six LBL objects selected. The black points are archival data while the respectively red and blue butterflies are the 2FGL and IFHL 
measurements. Gray points are the AP Librae data presented in this work. 
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Table Bl. Proposed LBL-type objects for VHE obser vations. The 2FGL nam e is given in the first column with the position in the 
second and third. The redshift measurement taken from jAckermann et al.l2Ql ih or jShaw et alj2Q13h is reported in the fourth column. 
The name (column 5) of the counterpart associated with the 2FGL source was found in the 2LAC catalog. The last column in the mane 
of the best suited instrument for observations. The sources ai‘e ranked by predicted flux above 200 GeV. 


2FGL name 


“J2000 

<5j2000 

redshift 

Association 

Instruments 

2FGL J1719.3+1744 

17h 

19’" 13.05” 

17° 45' 06.4" 

0.137 

PKS 1717+177 

VERITAS/MAGIC 

2FGL J0617.6-1716 

06'’ 

17’" 33.67” 

-17° 15' 22.8" 

0.098 

CRATES J061733.67-171522.8 

H.E.S.S. 

2FGL J0738.0+1742 

07*’ 

38’" 07.39” 

17° 42' 19.0" 

0.424 

PKS 0735+17 

VERITAS/MAGIC - H.E.S.S. 

2FGL J1559.0+5627 

15*’ 

58’" 48.29” 

56° 25' 14.1" 

0.3 

TXS 1557+565 

VERITAS/MAGIC 

2FGLJ1150.1+2419 

11 *’ 

50”" 19.21” 

24° 17' 53.8" 

0.2 

B2 1147+24 

VERITAS/MAGIC 

2FGL J0712.9+5032 

07*’ 

12’" 43.68” 

50° 33' 22.7" 

0.502 

GB6 J0712+5033 

VERITAS/MAGIC 









